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Abstract. Two different approaches are formulated to analyze two-dimensional quantum 
models which are not amenable to standard separation of variables. Both methods are essen- 
tially based on supersymmetrical second order intertwining relations and shape invariance - 
two main ingredients of the supersymmetrical quantum mechanics. The first method ex- 
plores the opportunity to separate variables in the supercharge, and it allows to find a part 
of spectrum of the Schrodingcr Hamiltonian. The second method works when the standard 
separation of variables procedure can be applied for one of the partner Hamiltonians. Then 
the spectrum and wave functions of the second partner can be found. Both methods are 
illustrated by the example of two-dimensional generalization of Morse potential for different 
• values of parameters. 
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\0 \ 1 Introduction 

^ , The exactly solvable models in quantum mechanics are of special interest during many years 

both by methodological and practical reasons. By now, the main achievements were related to 
one-dimensional Shrodinger equation. Indeed, a list of exactly solvable one-dimensional prob- 
lems (Harmonic oscillator, Coulomb, Morse, PoschRTeller potentials etc.) was obtained by 
an algebraic procedure in the framework of factorization method [T] in the middle of last cen- 
tury. This method was reproduced rather recently in supersymmetrical quantum mechanics 
approach [2] initiated by the seminal papers of E. Witten [3]. More of that, this approach gave 
many new exactly solvable potentials which were obtained as superpartners of "old" exactly 
solvable models. It is necessary to mention also the important paper of L. Gendenstein [4], 
where the new fruitful notion of shape invariance was introduced. For the sake of truth, more 
than a century ago the so called Darboux transformation [5] for Sturm-Liouville equation was 
well known among mathematicians. Its application to a specific Schrodinger-like equation is 
actually equivalent [6[ |7j to the factorization method. 

The situation is much worse for two-dimensional quantum mechanics. The only regular 
method to solve analytically the Schrodinger equation is well known method of separation of 
variables [8]. This method replaces the two-dimensional problem by a pair of one-dimensional 
problems. It can be used for very restrictive class of models. Full classification of models 
which allowed separation of variables was given by L.P. Eisenhart [9j: four possibilities exist - 
Cartesian, polar, elliptic and parabolic coordinates. The general form of potentials amenable to 
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separation of variables is known explicitly up to arbitrary functions of one variable. And analyti- 
cal solution is possible only if these functions belong to the list of exactly solvable potentials. All 
these Hamiltonians H are integrable: the symmetry operator R of second order in derivatives (in 
momenta) exists: [H, R] = 0. Besides models amenable to separation of variables, the class of so 
called Calogero-like models [10] is known as well. They describe the specific forms of pairwise 
interaction of N particles on a line, and they are solvable by means of special transformation of 
variables which leads to a separation of variables. The problem is in analogous state for higher 
dimensions of space. 

Intermediate class of models - quasi-exactly-solvable (QES) potentials (or, partially sol- 
vable) - became interesting during last years. This notion concerns models for which only a part 
of spectrum and corresponding wave functions can be found analytically. In one-dimensional 
quantum mechanics a lot of such models were built with some hidden algebraic structure 
The supersymmetrical approach also gave some new QES potentials [12| . 

Thus, the search of new approaches to solution of nontrivial two-dimensional quantum me- 
chanical models seems to be of current importance. It was already mentioned above that super- 
symmetrical quantum mechanics provides both new ways to derive some old results and inter- 
esting method to obtain new ones. In this paper we shall focus on the case of two-dimensional 
Schrodinger equation. Namely, we shall present two procedures of using of the supersymmet- 
rical intertwining relations with supercharges of second order in derivatives as procedures of 
SUSY-separation of variables. 

In Section [2] the general form of two-dimensional supersymmetrical quantum mechanics with 
second order supercharges will be formulated. Section [3] presents the first procedure of SUSY 
separation of variables where variables are separated in the supercharge. It leads to QES models, 
and the specific model of two-dimensional Morse potential illustrates this method. In Section U] 
the second procedure of SUSY separation of variables is given where variables are separated in 
one of partner Hamiltonians. In the case of the same Morse model, but with particular values of 
parameter, it allows to solve the model completely, i.e. to find analytically the whole spectrum 
and all wave functions. 



2 Two-dimensional SUSY quantum mechanics 

Direct generalization of one-dimensional Witten's SUSY quantum mechanics to the arbitrary 
dimensionality d of space was formulated in [61 [13]. The Superhamiltonian included (d + 1) 
matrix components of different matrix dimensionality, and these components are intertwined by 
components of supercharge - operators linear in derivatives. In particular, in the case of d = 2 
two scalar Hamiltonians and a 2 x 2 matrix Hamiltonian are intertwined 

# (o) Hlk ^ m 
it pf 

where 

q f = Tdl + (dtWm Pt = ZikqJ, H(°\+ = 4 H ki' H {0) Pt = P + k H H- 

Some physical problems were considered in this framework. For example, the spectrum of the 
Pauli operator describing spin 1/2 fermion in the external electrostatic and magnetic field was 
investigated |14 4 115 1 fT6| . However, the following natural question arises: is it possible to avoid 
matrix Hamiltonians from the scheme? Any attempt to intertwine two scalar Hamiltonians 
by means of first order operators leads to potentials with standard separation of variables [T7] 
which are not interesting for us here. 
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The nontrivial way to avoid matrix Hamiltonians lies in the framework of polynomial SUSY. 
The latter was used for the first time [7] in one-dimensional SUSY quantum mechanics. 
In two-dimensional context it was proposed in |19j . where a pair of scalar two-dimensional 
Hamiltonians H^°\ iJW was intertwined by second order operators 

h (o)q+ = Q+H (i) j q-hW = hWq~, 

where the Hamiltonians have the Schrodinger form 

H W = -d 2 + V^ 1 \x). 

As for the intertwining supercharges Q^, the first naive idea is to choose reducible (factorized) 
supercharge Q + = qfqT 

it Pf 

<=> H lk <=> m 

W Hlk <=> HW 
of Pf 



It is too naive |19j . since this construction leads to Hamiltonians amenable to separation of 
variables in polar coordinates 

V(x) = a 2 p 2 + ±F(<p). 

P 2 

The second idea is to choose Q + = qfUi^qZ with some unitary twist by constant matrix Ui k 

if Pf 
H®4=> H ik <^F(°) 



tfW <=> U u H lm ul k <=> tf« 

?i Pi 

Some QES models were obtained by this trick |20j . 
The most general form of second order supercharges 



Q + = 9ik(md k + Ci{x)di + B(x), Q- = (Q+)t 

leads to a complicate system of nonlinear second order differential equations for functions gi k , 
Ci, B, and potentials V^ 0,1 \x). Its general solution is impossible, but some particular solutions 
were found [191 EI]- The simplest choice g% k {x) = 5i k gives the separation of variables in polar 
coordinates. The Lorentz form gi k {x) = diag(l,— 1) does not lead to separation of variables, 
and some particular solutions were found [19] . 

Here we focus on gi k {x) = diag(l, —1). In this case, the system is simplified essentially. New 
variables x± = %\ ±x 2 are useful together with x\, x 2 . Using the intertwining relations, one can 
prove that new functions C± depend on one variable only 

C + = C 1 -C 2 = C+(x+), C^ = C 1 +C 2 = C-{xJ), x ± = Xl ± x 2 . 

The general solution for Lorentz metric can be provided by solving the only equation 



d.(C.F) = -d+(C+F), 
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where new useful function is F = Fi(x + +x-)+F 2 (x + — x-). Thus, the equation is the functional 
differential equation, and no regular procedure of its solution is known. 

The required potentials V^ 0,1 \x) and the function B(x) are expressed in terms of C± and F\^ 2 

= ±\(C' + + C'_) + \{C 2 + + Cl) + \ (F 2 (x + - x.) - F 1 {x + + x.)) , 

B = j (C+C_ + + x-) + F 2 (x + - x_)) . 

A variety of such pairs of potentials was found in [19] . 

3 SUSY-separation of variables I: QES models 

The first variant of SUSY-separation of variables is realized when the Hamiltonian H does not 
allow standard separation of variables, but the supercharge Q + does allow [22} I23j. The general 
scheme is the following. Let's suppose that we know zero modes of Q + 

Q+n n (x) = 0, n = 0,l,...,N, Q + n{x) = 0. 

The intertwining relation 

ff(°)Q+ = Q+hM 

obey the important property: the space of zero modes is closed under the action of : 

H {1) n(x) = cn(x). 

If the matrix C is known, and if it can be diagonalized 

BC = AB, A = diag(A ,Ai,...,A7v), 
the eigenvalues of can be found algebraically 

H^(M{x)) = A(M(x)). 
Thus, for realization of this scheme we need 

- to find zero modes Cl n (x); 

- to find constant matrix B, such that BC = KB. 

As for zero modes, they can be obtained by using the special similarity transformation (not 
unitary!), which removes the terms linear in derivatives from Q + 

q+ = e -x(-)Q+ e +x(-) = 32 _ d 2 + I( Fl ( 2xi ) + F 2 (2x 2 )), 
X{x) = ~(J C + (x+)dx + + j C^(x^)dx^ 

Now, q + allows separation of variables for arbitrary solution of intertwining relations, and we 
obtain the first variant of new procedure - SUSY-separation of variables. Similarly to the 
conventional separation of variables, separation of variables in the operator q + itself does not 
guarantee solvability of the problem. 

The next task is to solve two one-dimensional problems 

-df - ^F 1 {2x 1 ) \ rj n (xi) = e n i] n (xi), 

-d\ + ^F 2 {2x 2 ))pn{x2 S j = e n p n (x 2 ). 
Three remarks are appropriate now. 
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Remark 1. The same similarity transformation of does not lead to operator amenable to 
separation of variables. 

Remark 2. The normalizability of £l n has to be studied attentively due to non-unitarity of the 
similarity transformation. 

Remark 3. We have no reasons to expect exact solvability of the model, but quasi-exact- 
solvability can be predicted. 

As for the matrix B, it must be found by some specific procedure. Such procedure was used 
in example which will be presented below. 

In principle, the first scheme of SUSY-separation of variables can be used for arbitrary models 
satisfying intertwining relations by supercharges with Lorentz metrics. The list of solutions of 
intertwining relations is already rather long, and it may increase in future. The main obstacle 
is analytical solvability of one-dimensional equations, obtained after separation of variables in 
the operator q + . 

Below we describe briefly such a model which can be considered as the generalized two- 
dimensional Morse potential 

C+ = 4aa, C- = Aaa coth ■ 



2 ' 

fiixi) = -M2 Xi ) = -A (e- 2ax > - 2e~^) , i = 1, 2, 



4 

'ax. 



■ctX2\ 



VK°)>« = a 2 a(2a T 1) sinh" 2 (— J + 4a 2 a 2 + A [e~ 2aXl - 2e~ aXl + e~ 2aX2 - 2e 

where A > 0, a > 0, a is real. 

To explain the name, we present the potential in the form 

V(x) = VMorscOl) + ^MorseO^) + v(xi, X 2 ), 

where first two terms are just one-dimensional Morse potentials, and the last term mixes vari- 
ables Xi, X2- 

The solutions of one-dimensional Schrodinger equations are well known [23], and the zero 
modes can be written 1221 [231 as 



a n (x) = (^Jit^j exp f-^y^) (^) Sn F(-n,2s n + l;^)F(-n,2 Sn + 1;£ 2 ), 

t 2^1 v 7 ! 1 

& = exp(-axi), s n = n-->0. 

a a 2 

The conditions of normalizability and of absence of the "fall to the center" are 



/ 1 1 \ VA 1 

aG h"I"Wi J' s n = — -n-->-2a>0. 

To obtain the matrix C explicitly, one must act by on £l n . The matrix turns out to be 
triangular, and therefore, the energy eigenvalues coincide with its diagonal elements 

Ek = c kk = -2(2aa 2 s k - e k ). 

To find a variety of wave functions is a more difficult task. For that it is necessary to find 
all elements of C and all elements of matrix B. The recurrent procedure for the case of two- 
dimensional Morse potential was given in |22[ [23] . This variety can be enlarged by means of 
shape invariance property [25] of the model 

tf (0) (£;a) = # (1) (x;a) + H(a), 5 = a -1/2, H(a) = a 2 (4a - 1). 
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Similarly to one-dimensional shape invariance, each wave function constructed by SUSY-sepa- 
ration of variables leads to a set of additional wave functions 



M — 1 



-t) 



~{a)Q a 



1 



M - 1 
2 



M - 1 



Analogous approach works for the two-dimensional generalization of Poschl-Teller model [20 
and for some two-dimensional periodic potentials [26] . 



4 SUSY-separation of variables II: exact solvability 

Among all known solutions of two-dimensional intertwining relations with second order super- 
charges a subclass exists |27| . where one of intertwined Hamiltonians is amenable to standard 
separation of variables due to specific choice of parameters of the model. Its superpartner still 
does not allow separation of variables. 

The scheme will be described below for the same specific model which is two-dimensional 
generalization of Morse potential 



VK0),(1) = a 2 fl(2a -p x) sinh -2 fO^-\ + 4a 2 a 2 + A [ e -2a Xl _ 2& -ax 1 + e ~2ax 2 _ 2& - 



ax2 



Let's choose ao = —1/2 in order to vanish the mixed term in . Then allows the 

conventional separation of variables. Moreover, after separation of variables each of obtained 
one-dimensional problems is exactly solvable. We met just this one-dimensional problem above 
in a different context. 

The discrete spectrum of this one-dimensional model is 



2 2 v A 1 

e n = -a s n , s n = n-->0, n = 0,1,2, 

a I 

Wave functions are expressed in terms of degenerate hypergeometric functions 

r) n (xi) = exp (&) Sni? (- n > 2s n + C» = ^^exp(-axi). 

Due to separation of variables, the two-dimensional problem with H^(x) is exactly solvable. 
Its energy eigenvalues are 

En,m — E mn — 6 n + £ m , 

being two-fold degenerate for n ^ m. The corresponding eigenfunctions can be chosen as 
symmetric or (for n ^ m) antisymmetric combinations 

*^^ A (^) = r} n {x 1 )-q m (x2) ±r] m (xi)r] n (x 2 ). 

Our aim here is to solve completely the problem for H^°\x) with ao = —1/2. The main 
tool is again the SUSY intertwining relations, i.e. isospectrality of and but up to zero 
modes and singular properties of Q^. In general, we may expect three kinds of levels of H^°\x): 
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(i) The levels, which coincide with E nm . Their wave functions can be obtained from by 
means of Q + . 

(ii) The levels, which were absent in the spectrum of (x), if some wave functions of (x) 
are simultaneously the zero modes of Q~ . Then the second intertwining relation would 
not give any partner state among bound states of H^(x). 

(iii) The levels, which were also absent in the spectrum of H^\x), if some wave functions 
of H(°\x) become nonnormalizable after action of operator Q~. 

We have to analyze these three classes of possible bound states of H^°> one after another. 

(i) The first SUSY intertwining relation gives the two-fold degenerate wave functions of H ^ 
with energies E nm : = Q + tyQ . But Q + includes singularity on the line x\ = x 2 , therefore 
the normalizability of m depends crucially on the behavior of m on the line £1 = £2- One 

can check that only antisymmetric functions survive, i.e. only symmetric survive. This 
fact can be demonstrated [2Tj both by direct calculation and by indirect method - by means of 
symmetry operator 

The indirect method explores that the symmetry operator R(°> = Q~Q + for ao = —1/2 can 
be written in terms of one-dimensional Morse Hamiltonians hi(x\), h 2 {x 2 ) 

= (hxixt) - h 2 {x 2 )f + 2a 2 (/i 1 (x 1 ) + h 2 {x 2 )) + a 4 . 

Therefore, 

RM^l A Jx) = r n>m *£ n Jx), r n , m = a%n - raf - l][(s n + s m f - l] , 

and 

||^g)5 ||2 = { ^0)A j + j g M )= rnm \\*P A || 2 . 

For n = m, wave functions vanish identically by trivial reasons. It is clear now that 

wave functions $2 also vanish. For all other n, m, functions $2 have positive and finite 
norm, and there is no degeneracy of these levels. 

(ii) These possible bound states of are the normalizable zero modes of Q~ . The variety 
of such zero modes is known from [22] : they exist only for positive values of a 

( l 1 

a G - H = , +00 

which does not contain the value ao = —1/2. Thus, no normalizable bound states of this class 
exist for H^°\ 

(iii) We have to study an opportunity that Q~ destroys normalizability of some eigenfunc- 
tions of H(°\ It could occur due to singular character of Q at x\ = x 2 . The analysis was 
performed [27] in suitable coordinates. It shows that Q~ is not able to transform normalizable 
wave function to nonnormalizable. Therefore, the third class of possible wave functions 
does not exist too. 

Summing up, the spectrum of with ao = —1/2 consists only of the bound states with 
energies E nm for \n— m\ > 1. This spectrum is bounded from above by the condition of positivity 
of s n , s m : n,m < \f~Aj a — 1/2. The corresponding wave functions are obtained analytically |27] . 

The results above can be expanded to the whole hierarchy of Morse potentials with = 
— {k + l)/2 with k = 0,1, . . . by means of shape invariance property. Let's denote elements of 



8 



M.V. Ioffe 



the hierarchy as H^°\x; a k ), H^\x;a k ). All these Hamiltonians are also exactly solvable due 
to shape invariance of the model 

(x; afc_i) = H^(x; a k ), k = 1,2, 

This means that the following chain (hierarchy) of Hamiltonians can be built 

H^\x;a )^H ( - \x;a ) = H^(x;a 1 )^H ( - \x;a 1 ) = ■ ■ ■ -=- H^\x; a fc _i) 

= H^{x;a k ) + H^{x;a k ), 

where the sign -f- denotes intertwining by Q ± (ai). 
In the general case, the functions 

(if normalizable) are the wave functions of (x; a k ) with energies E n m = — a 2 (s 2 + s^J. The 
symmetries of wave functions alternate and depend on the length of chain. This is true but up 
to zero modes of operators Q + . 

It is necessary to keep under the control normalizability of VP and zero modes of Q + . This 
control is performed algebraically by means of identity, which must be fulfilled up to a function 
of H 

R^(a k ) = R^\a k _ 1 ). 

Actually, the following equation holds: 

Q~(a k )Q + (a k ) = Q+(a fc _ 1 )Q-(a fc _ 1 ) + a 2 (2k + 1) [2H^(x; a k -i) + a 2 (2k 2 + 2k + 1)] . 

These relations allowed to evaluate the norms of wave functions. The result is the following. 
The spectra of Hamiltonians H^°'(x;a k ) are not degenerate. They consist of the bound states 
with energies -E n ,m 5 with indices |n — m\ > k + 2, and their wave functions (x; Ofc) were 

given analytically above. 
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